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FOREWORD  ' 

The  work  reported  here  was  performed  in  the  Mechanics  and  Surface 
Interactions  Branch,  Nonmetalllc  Materials  Division,  Air  Force  Materials 
Laboratory,  Wright-Patterson  AFB,  Ohio.  H.  Thomas  Hahn,  AFML/MBM,  was 
the  principal  Investigator  and  R.  Y.  Kim  of  the  University  of  Dayton 
Research  Institute  was  responsible  for  the  experimental  phase  of  the 
project. 

The  authors  wish  to  acknowledge  R.  Esterline  and  R.  Cornwell  for  the 
preparation  and  testing  of  specimens.  Sincere  appreciation  Is  extended 
to  S.  Dastin  of  the  Grumman  Aerospace  Corporation  for  supplying  some  of 
the  specimens  tested. 
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SECTION  I 
INTRODUCTION 

Susceptibility  of  polymers  to  swelling  due  to  moisture  absorption 
leads  to  the  presence  of  residual  stresses  in  the  composite  laminates 
which  utilize  such  materials  as  matrix.  In  composite  laminates  a non- 
uniform  moisture  distribution  is  not  the  only  cause  of  the  swelling 
stress;  even  under  a uniform  distribution  the  residual  stresses  result 
from  the  incompatible  swelling  behavior  between  the  matrix  and  fibers. 

The  swelling  stresses  thus  induced  may  lead  to  the  creation  of  micro- 
cracks especially  under  transient  conditions  and  further  degrade  the 
strength  properties  of  composites.  At  the  same  time,  however,  these 
stresses  generally  negate  the  curing  stresses  and,  hence,  the  combined 
effects  may  be  beneficial.  Thus,  it  is  important  to  know  not  only  the 
moisture  diffusion  behavior  but  also  the  swelling  behavior,  in  order  to 
properly  characterize  the  material  property  change  resulting  from  the 
moisture  absorption. 

Although  the  interfacial  effects  are  not  negligible,  many  of  the 
moisture-related  characteristics  of  composites  are  known  to  be  traceable 
back  to  those  of  the  matrix  material.  For  example,  moisture  diffusion 
in  graphite/epoxy  composites  is  described  by  the  Fick  equation  with  a 
diffusion  coefficient  which  can  be  predicted  from  the  matrix  properties 
(Reference  1).  The  matrix/ interface-control  led  properties  degrade 
substantially  at  the  high  temperature  wet  conditions  (Reference  2)  and 
the  glass  transition  temperature  of  composite  is  lowered  similarly  to 
that  of  matrix  (Reference  3).  The  equilibrium  swelling  strains  of 
composite  and  of  matrix  can  be  related  to  each  other  through  a proper 
analysis  (Reference  4). 

In  this  report  some  experimental  results  for  a graphite/epoxy  com- 
posite (AS/3501-5)  are  presented  to  support  the  validity  of  the  swelling 
strain  model  proposed  in  Reference  5.  The  model  is  basically  derived 
from  the  result  of  a micromechanics  analysis,  but  it  can  be  modified 
slightly  to  accomodate  the  experimental  observations.  Both  lamina  and 
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laminate  swelling  strains  are  analyzed,  and  residual  stresses  calculated, 
by  this  model.  A difference  between  the  absorption  and  desorption  rates 
is  explained  in  terms  of  the  residual  stresses  and  stress  dependency  of 
diffusion.  In  addition,  warpage  of  unsymmetric  laminates  is  analyzed 
and  the  matrix  cracking  is  discussed.  Finally  presented  is  a relation- 
ship between  the  temperature  and  the  relative  humidity  which  results  in 
a residual  stress-free  state  for  the  graphite/ epoxy  composite. 
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SECTION  II 

THEORETICAL  BACKGROUND 

1.  HYGROTHERMAL  STRAINS  OF  UNIDIRECTIONAL  LAMINAE 


It  has  been  shown  in  Reference  5 that  the  swelling  strains  of  a 
composite  reinforced  with  moisture-insensitive  fibers  are  given  by 


1 + V 

— r s(c  - c )H(c  - c ) 

3 V V 


(1) 


Here  is  the  Poisson's  ratio  of  matrix,  s and  c are  the  specific 
gravity  and  the  moisture  content,  respectively,  of  the  composite,  and 
Cy  represents  the  amount  of  moisture  entrapped  in  the  voids 


c = V /s 
V V 

where  v^  is  the  void  content.  H(*)  is  the  Heaviside  step  function. 


(2) 


In  reality,  the  volume  additivity,  which  has  been  assumed  in  the 
fore-going  derivation,  may  not  hold,  but  the  swelling  strain  is  still 
linearly  related  to  the  moisture  concentration  above  a certain  thresh- 
hold  level.  In  this  case  Equation  1 can  be  generalized  to 

e”  = (c  - c^)H(c  - c^)  (3) 

U 

Equation  3 is  graphically  shown  in  Figure  1.  Note  that  is  the  trans- 
verse swelling  coefficient  to  be  determined  experimentally  and  that  c^ 
can  now  be  regarded  as  a threshhold  value  of  c below  which  no  appre- 
ciable swelling  occurs. 


In  the  presence  of  temperature  changes,  the  resultant  nonmechani 
cal  strains  of  a unidirectional  composite  are  given  by 


e 


L 


(4) 
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If  thermal  expansion  is  linear  with  temperature  change,  then  the  thermal 
strains  are 

-I  - <*’ 

In  the  subsequent  analyses,  the  material  behavior  is  assumed  elastic 
and  the  swelling  strains  independent  of  temperature.  For  convenience, 
the  thermal  strains  are  measured  in  the  dry  state  as  the  temperature  is 
changed  from  the  initial  stress-free  reference  to  the  final  temperature 
and  then  the  swelling  strains  are  superimposed  at  the  final  temperature. 
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2.  HYGROTHERMAL  STRAINS  OF  MULTIDIRECTIONAL  LAMINATES 

The  laminate  constitutive  relations  within  the  framework  of  the 
classical  laminated  plate  theory  are  given  In  Reference  7.  The  salient 
features  of  the  theory  are  briefly  described  here. 

The  constitutive  relations  of  a constituent  ply  In  a laminate  can 
be  written  as 


r.  = C..(e.  - e.)  + C (e  -el 

^ J j iA'  A ®a' 


^A  = - ®b) 


(6) 

(7) 


where  are  the  total  strains. 


The  subscripts  1 through  n are  associated  with  the  In-plane  coor- 
dinates and  the  subscripts  A through  C with  the  out-of-plane  coordinates. 
C's  are  the  stiffnesses  at  the  final  state  of  Interest. 

According  to  the  classical  laminated  plate  theory,  are  assumed 
to  be  constant  through  the  thickness  and  the  total  strains  are  given 
by 


(8) 


where  k^  are  the  curvatures  and  z Is  the  thickness  coordinate.  Therefore, 


Introducing  the  reduced  stiffnesses  Q 


1j 


-] 


Q = C - C .C  ' C 
ij  ij  lA  AB  Bj 

one  can  derive  the  laminate  constitutive  relations  as 


N.  = 

A €° 

+ B.  .k.  • 

■N^ 

+ T., 

1 

ij  J 

J 

1 

i] 

M,  = 

B,.e? 

+ D..k.  - 

-M^ 

+ R. 

i 

J 

U J 

1 

1 

®A  ^AB  ""b  ■ ^AB^Bj  ®j  ■ ^AB^Bj®  ''j 


(9) 

(10) 

(11) 

(12) 


^ ^ABSj®J  ^ ®A 
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N N 

Here  B^.j,  and  are  the  usual  laminate  stiffnesses,  and  and  M” 

are  the  nonmechanical  forces  and  moments,  respectively.  An  over  bar 
stands  for  the  average  through  the  thickness  h. 


( ) 


1 

si.  ‘ 


-h/2 


(13) 


and 


^h/2 


^iA^AB 


[ 1 , z]dz 


(14) 


The  nonmechanical  laminate  strains  and  curvatures  are  obtained 
by  putting  = Og  = 0 in  Equations  10  through  12: 


where 


= f7/(N^  - B D'^  M^) 

1 j IK  kn  n 

(15) 

,N  r^-1  „ oN 

N ' °ij  '“j  - ' 

(16) 

®A  = *A  + ^AR^R  ®-  • ^Ir^R- 

A A ABBjj  ABBj  j 

AB  Bj  j 

(17) 

F..  = A..  - B.,  D,'^B  . 

ij  ij  ik  kn  nj 

(18) 

The  in-plane  residual  stresses  in  each  ply  are  then  obtained  from 
Equation  6 with  the  aid  of  Equation  9: 


^ , oN  ^ . N 
Q.A^t  + k,  z - 
ij  J j 


(19) 


Moisture  diffusion  in  composite  laminates  has  been  shown  to  be 
described  fairly  well  by  the  Pick  equation  (Reference  1).  The  solution 
to  this  equation  for  a thin  laminate  undergoing  moisture  absorption  or 
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desorption  from  an  equilibrium  state  is  well  known  (for  example,  see 
Reference  8).  The  expressions  for  the  moisture  concentration  c and 
the  average  moisture  content  through  the  thickness,  c",  are  repeated  here 
for  convenience: 


c - c. 


V 


- i: 


j=0 


(-  1)^ 
(2j  + 1) 


cos 


(2i  + z 


exp  I- 


12 


2 2 \ 
1 -f  n Dt  \ 

/ 


(20) 


00  0 


00 


= 1-^  L 


exp 


(2i  + 1)^  n^Dt 


(21) 


Here  D is  the  diffusivity,  t the  time,  and  the  subscripts  0 and  ” stand 
for  the  initial  and  the  final  equilibrium  state,  respectively.  Note  that 


c„  = c and  c. 


'O' 


The  through-the-thickness  distribution  of  the  unconstrained  swelling 
strain  follows  upon  substitution  of  Equation  20  into  Equation  1 or  3,  in 
general.  The  final  laminate  strains  and  curvatures  due  to  swelling  are 
then  obtained  by  solving  Equations  15  through  18. 


H 

For  symmetric  laminates  the  nonmechanical  moments  vanish  because 
c is  also  symnetrically  distributed  with  respect  to  the  mid-plane. 
Furthermore,  since  c depends  on  z/h  and  also  since  the  nondimensional 
time  Dt/h  can  be  expressed  in  terms  of  c from  Equation  21,  the  laminate 
swelling  strains  can  be  shown  to  depend  only  on  the  average  moisture 
content  independently  of  the  thickness  and  diffusivity.  Thus,  all  lami- 


nates of  the  layup 


^ "2  "3 


], 


will  follow  the  same 


relationship  between  the  swelling  strains,  both  in-plane  and  out-of-plane, 
and  the  average  moisture  content  c,  as  long  as  the  ply  ratio  0^:02:02.’  — 
remains  the  same.  Since  the  relationship  is  independent  of  the  diffusivity, 
so  it  is  of  the  temperature  insofar  as  the  equilibrium  moisture  content 
remains  the  same.  The  aforementioned  result  is  useful  in  that  the 
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swelling  strains  of  symmetric  laminates  can  be  determined  within  the 
shortest  time  by  testing  the  thinnest  possible  laminate  at  the  highest 
allowable  temperature. 


For  unidirectional  laminae.  Equation  15  for  the  swelling  strain 
reduces  to 


or 


H 


h 


-h/2 


(c  - c )H(c  - c )dz 

V V 


H 

Qj 


V C 


(22) 


(23) 


where  y is  the  fraction  of  the  moisture  content  that  contributes  to 
swelling.  If  the  moisture  profile  through  the  thickness  is  as  shown 
in  Figure  2,  then  y is  the  shaded  area  divided  by  the  total  area  under 
the  moisture  profile. 

In  an  absorption  test  it  is  possible  that  c>Cy  everywhere.  In  this 
case  Equation  23  reduces  to 


oH  H ““  ““ 

*T  ^ (c  - c^)  if  c > (24) 

where  c^  is  the  value  of  c when  c reaches  c^  everywhere.  Thus  the 
swelling  coefficient  and  the  threshhold  value  c„  can  be  determined 
from  the  linear  region  of  the  plot  of  e-f  versus  c.  However,  the 
same  is  not  true  for  the  desorption  test.  The  relationship  between 
and  Cy  is  shown  in  Figure  3. 
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THICKNESS 


Figure  2.  Fraction  of  Moisture  Content  Which  Contributes  to  Swelling 
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SECTION  III 

EXPERIMENTAL  CORRELATION 
1.  EXPERIMENTAL  PROCEDURE 

Two  types  of  experiments  involved  in  this  study  are  measurments 
of  swelling  strains  of  symmetric  laminates  and  of  warpage  of  unsymmetric 
laminates  as  a function  of  moisture.  The  material  system  chosen  is 
graphite/epoxy  AS/3501-5  except  as  indicated  in  Table  1. 

All  specimens  but  [0^/90^].|.  plates  were  cut  out  of  larger  panels 
which  were  fabricated  from  prepreg  tape  in  an  autoclave  according  to 
the  manufacturer's  recommended  cure  cycle.  Layups  and  dimensions  of 
the  specimens  are  listed  in  Table  1.  Limited  variation  of  dimension 
was  allowed  to  study  any  possible  size  effect.  All  specimens  were  stored 
in  a vacuum  oven  at  71°C  for  at  least  one  week  prior  to  the  test  in 
order  to  dry  out  any  pre-existing  moisture  absorbed  during  the  period 
of  storage. 

Thermal  strains  were  measured  by  strain  gages  (Micro-Measurements 
CEA-00~250  UW-120)  in  accordance  with  the  fused  silica  compensating 
technique.  The  specimen  temperature  was  monitored  with  a type  K 
thermocouple  directly  mounted  on  the  specimen  surface.  The  specimen 
heating  was  achieved  in  a convection  type  oven  at  the  rate  of  2.8°C/ 
min.  Strain  readings  were  taken  over  two  consecutive  temperature 
cycles  ranging  from  23  to  182°C. 

Moisture  conditioning  was  achieved  at  82°C  in  distilled  water 
for  absorption  and  in  vacuum  for  desorption.  The  time  interval  be- 
tween the  first  and  second  conditioning  was  about  six  months  and  the 
data  were  analyzed  independently  for  each  cycle. 


TABLE  1 
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Swelling  strains  were  measured  by  either  a micrometer  of  gage  length 
2.54  cm  or  a caliper  of  gage  length  20.3  cm.  The  resolutions  of  the 
micrometer  and  the  caliper  are  0.00254  and  0.0127  mm,  respectively. 

Strain  measurements  were  taken  at  four  locations  for  the  thickness 
expansion  and  at  three  locations  for  the  transverse  and  axial  directions. 
The  longitudinal  strain  of  unidirectional  lamina  and  the  in-plane 
strains  of  multidirectional  laminates  are  not  reported,  since  they  are 
negligibly  small . 

The  unsymmetric  plates,  supplied  by  the  Grumman  Aerospace  Corp., 
we’^a  fabricated  under  two  different  cure  processes:  a cooldown  rate  of 
0. 55°C/min  to  121°C  followed  by  a post  cure  of  3 hrs  at  187.8  ± 5.6°C 
for  the  specimens  US-1  through  7,  and  a cooldown  rate  of  2.8°C/min 
followed  by  a post  cure  of  8 hrs  at  176.7  ± 5.6°C  for  the  specimens  UF- 
1 through  7. 

The  plates  were  initially  warped  due  to  the  residual  curing  stresses. 
The  deflection,  defined  as  the  relative  displacement  of  the  center  of 
the  plate  when  placed  on  a flat  surface,  was  monitored  as  a function 
of  moisture  content  by  a travelling  telescope  with  resolution  of  0.01 
mm.  The  same  moisture  conditioning  as  in  the  case  of  in-plane  strain 
measurements  was  employed  and  the  tests  were  terminated  when  the 
deflection  nearly  vanished. 


2.  RESULTS  AND  DISCUSSION 


The  transverse  swelling  strains  obtained  from  the  unidirectional 

specimens  are  shown  in  Figure  4.  During  absorption  both  specimens 

exhibit  almost  the  same  behavior  independently  of  the  test  sequence. 

The  threshhold  value  c^  above  which  the  swelling  strain  increases 

linearly  with  the  moisture  content  is  seen  to  be  approximately  0.9 

percent.  The  slope  of  this  linear  portion  is  the  swelling  coefficient 
H 

oj  and  the  intersection  of  its  extension  with  the  abscissa  yields  c^: 


Oj  = 0.59  and  c^  = 0.4% 
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Figure  4.  Transverse  Swelling  Strain  of  [0]gj  Laminate 
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; It  is  interesting  to  note  that  the  measured  swelling  coefficient  is  lower 

than  the  theoretical  value  of  0,72  that  follows  from  Equation  1 and  the 
I material  properties 

I 

I 

I V = 0.35  and  s = 1.6 

' m 

u 

Once  oij  and  c^  are  known,  the  swelling  strain  can  be  calculated 
over  the  entire  range  of  the  moisture  content.  The  solid  curve  in  Figure 
1 4 is  the  result  of  this  calculation  when  using  the  equilibrium  moisture 

content  of  1.7%. 

i During  desorption,  however,  the  analytical  result  does  not  agree 

I well  with  the  experimental  data.  Only  the  specimen  0-2  behaves  close 

to  the  prediction.  The  other  specimen  follows  rather  a linear  relation- 
ship than  a bilinear  one.  The  same  behavior  is  repeated  in  the  second 
desorption  test. 

A difference  between  the  absorption  and  the  desorption  behavior 
appears  also  in  the  rate  of  change  of  the  moisture  content.  Figure  5. 

Four  distinctive  features  seem  worthy  of  elaboration:  (1)  The  equili- 
brium weight  change  increases  after  each  moisture  conditioning;  (2) 

There  is  a weight  loss  and  contraction,  albeit  small,  after  desorption; 

(3)  The  initial  absorption  of  virgin  specimens  facilitate  the  subsequent 
diffusion;  (4)  The  desorption  takes  place  at  faster  rate  than  the  absorp- 
tion. Furthermore,  the  linear  portion  of  the  weight  change  versus  (time) 

1/2 

' relationship  is  longer  in  absorption  than  in  desorption. 

The  first  two  observations  indicate  the  possibilities  of  microvoid 
formation  and  of  leaching  out  of  low  molecular  weight  material  (Reference 
9).  The  third  points  toward  the  additional  possibility  of  microcracking 
during  the  first  cycle  which  may  open  up  the  passageway  for  the  subsequent 
moisture  diffusion  but  do  not  necessarily  entrap  moisture.  The  last 
observation  as  well  as  the  poor  experimental  correlation  in  Figure  4 can 
be  explained  in  terms  of  the  residual  stresses  due  to  swelling  as  follows. 
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Calculation  of  the  residual  stresses  requires  the  elastic  properties 
to  be  known  at  the  condition  of  interest.  Since  the  moisture  concentra- 
tion varies  through  the  thickness  during  the  transient  stage  of  the 
moisture  diffusion,  the  elastic  properties  must  be  available  for  the 
entire  range  of  the  moisture  concentration.  Table  2 lists  the  available 
elastic  properties  at  the  room  temperature  dry  (RTD)  and  room  temperature 
wet  (RTW)  conditions.  Since  difference  between  the  two  sets  of  properties 
is  rather  small,  however,  the  average  properties,  as  listed  in  Table  3, 
are  used  in  the  analysis.  Table  3 also  lists  the  thermal  expansion 
coefficients  measured  in  the  present  study.  It  is  noted  in  passing  that 
the  thermal  expansion  behavior  was  reproducible  over  the  two  complete 
cycles  between  the  room  temperature  and  182°C. 

The  transverse  residual  stress  at  the  lower  surface  of  each  ply  is 
calculated  during  absorption  and  shown  in  Figure  6(a).  The  numbers 
represent  the  ordinal  position  of  each  ply  from  the  outside  inward. 

Thus,  the  curve  numbered  1 is  for  the  outside  surface.  The  stress  changes 
from  compression  to  tension  as  one  goes  inward  from  the  outside  surface. 

The  compression  near  the  boundary  slows  down  the  diffusion  and  the  more 
realistic  moisture  profile  will  be  as  shown  in  Figure  6(b)  (References 
8 and  10).  At  the  same  moisture  content,  therefore,  the  actual  swelling 
strain  will  be  smaller  than  predicted,  the  difference  being  proportional 
to  the  net  change  in  the  shaded  area.  Note  that  this  difference  appears 
only  up  to  when  the  moisture  concentration  at  the  mid  surface  is  less 
than  the  threshhold  value. 

During  desorption,  on  the  other  hand,  the  outside  ply  is  in  tension 
(Figure  7(a),  resulting  in  accelerated  desorption.  The  corresponding 
moisture  profile  in  Figure  7(b)  indicates  that  the  actual  swelling 
strain  increases  from  the  predicted  in  proportion  to  the  net  change  in 
the  shaded  area. 

Although  the  transverse  stress  at  the  outside  surface  exceeds  the  trans- 
verse strength,  no  visible  cracks  were  detected  by  microscopic  examinations 
of  up  to  x300  magnification.  The  reason  may  be  because  only  a small  region 
is  subjected  to  the  high  tensile  stress  while  the  rest  are  in  compression. 
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Figure  6(b).  Predicted  Moisture  Profile  (Solid  Curve)  and  Actual 
Moisture  Profile  (Broken  Curve)  During  Absorption 
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Figure  7(b).  Predicted  Moisture  Profile  (Solid  Curve)  and  Actual 
Moisture  Profile  (Broken  Curve)  During  Desorption 
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The  thickness  strains  of  the  [0]gj  and  [0/90^]^  laminates  are  shown 
in  Figures  8(a)  and  8(b),  respectively.  For  the  [0]gj  laminate,  the 
analytical  result  is  the  same  as  that  for  the  transverse  strain.  Large 
scatter  in  the  data  is  the  result  of  the  specimens  being  too  thin;  the 
thickness  is  in  the  range  of  1 mm.  An  error  in  reading  the  micrometer  can 
add  or  subtract  as  high  a strain  as  0.2  percent.  Thus,  the  data  should 
be  taken  as  an  indication  of  a trend  rather  than  as  the  material  property 
itself. 

The  absorption  and  desorption  behavior  of  the  [0/90g]^  laminate  is 
very  similar  to  that  of  the  [0]3i-  laminate.  Figure  9.  The  initially  ten- 
sile, transverse  residual  stresses  in  the  inside  plies,  which  result  from 
the  cure,  become  compressive  as  the  moisture  concentration  increases  in 
these  plies.  Figure  10(a).  In  Figure  10(a)  the  letters  L and  U stand  for 
the  lower  and  the  upper  surface,  respectively.  During  desorption  the 
outside  surface  is  subject  to  a tensile  transverse  stress  which  remains 
comparable  to  the  curing  stress  independently  of  the  moisture  loss.  Figure 
10(b). 


The  warpage  of  the  square  [0^/90^]^  laminate  is  calculated  from 
k!J  by  using  the  equation 


where  a is  the  length  of  one  side.  Both  analytical  and  experimental 
results  are  shown  in  Figure  11.  The  warpage  is  not  exactly  symmetric 
on  both  sides  and  only  the  maximum  values,  except  at  the  dry  state,  are 
given  for  the  reasons  described  later.  The  initial  warpage  is  due  to 
the  curing  stress  and  its  prediction  is  based  on  the  thermoelastic 
properties  in  Table  3. 


The  data  in  Figure  11  indicates  that  the  specimen  US-5  exhibits 
more  warpage  than  the  UF-5.  To  check  if  this  is  generally  representative 
of  the  difference  between  US  and  UF  specimens,  warpage  measurements 
were  taken  on  other  specimens  and  the  results  are  listed  in  Table  4. 
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Specimens  were  kept  in  water  at  82°C  for  one  year. 
Specimens  were  used  in  the  swelling  tests. 
Deflections  were  measured  in  the  dry  state. 

No  test  was  done. 


Figure  8(a).  Thickness  Swelling  Strain  of  Laminate 
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A11  specimens  except  those  numbered  4 and  6 had  been  kept  In  the  room 
environment  for  more  than  a year  and  are  believed  to  be  In  an  equilibrium 
state.  In  the  table  the  ply  angle  Is  measured  from  the  fiber  direction 
of  the  ply  that  was  In  contact  with  the  support  platen  during  cure.  The 
average  maximum  warpage  Is  1.7  nin  for  the  US  specimens,  compared  with  1.3 
mm  for  the  UF  specimens,  Indicating  the  validity  of  the  aforementioned 
observation. 

The  measured  warpage  Is  seen  to  be  greater  than  the  predicted. 
However,  the  slopes  are  almost  the  same.  Indicating  that  the  discrepancy 
Is  due  rather  to  the  cure  than  to  the  swelling. 

To  determine  the  cause  of  the  warpage  being  unsymmetric,  each 
specimen  was  polished  on  two  perpendicular,  adjacent  edges  and  trans- 
verse cracks  were  located  and  counted.  Examination  of  a few  specimens 
on  two  parallel,  opposite  edges  revealed  that  most  cracks  found  were 
through  the  entire  specimen  width.  A typical  crack  Is  shown  In  Figure 
12. 


Comparison  of  the  warpage  with  the  number  of  cracks  In  Table  4 
Indicates  that  the  less  contraction  Is  due  to  the  presence  of  cracks 
except  for  the  specimen  UF-1 . The  warpage  resulting  from  the  contraction 
of  the  plies  with  fewer  cracks  Is  larger  for  the  US  specimens  than  for 
the  UF  specimens.  This  Is  In  line  with  the  observation  that  the  average 
number  of  cracks  1n  the  US  specimens  (10)  Is  higher  than  that  In  the 
UF  specimens  (7)  since  both  conditions,  more  warpage  and  more  cracks, 
are  Indicative  of  higher  residual  stress.  For  some  reasons  unknown  the 
US  specimens  Invariably  have  more  cracks  In  the  90-deg  plies,  whereas 
the  UF  specimens  follow  a rather  random  trend  slightly  In  favor  of  the 
0-deg  plies. 

The  transverse  stresses  at  the  outside  and  mid  surfaces  are  plotted 
against  the  moisture  gain  In  Figure  13.  The  curing  stress  at  the  mid 
surface  Is  higher  than  that  at  the  outside  surface,  but  It  Is  only 
slightly  higher  than  half  the  RTD  strength. 
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Examination  of  Table  4 reveals  that  the  various  hygrothermal  condi- 
tioning has  little  to  do  with  the  number  of  cracks.  Figure  14  shows  the 
transverse  residual  stresses  at  the  room  temperature  as  the  equilibrium 
moisture  content  varies.  The  corresponding  warpage  is  seen  to  agree  well 
with  the  prediction  in  Figure  15.  In  the  room  environment  the  equilibrium 
moisture  content  was  found  to  be  about  0.6  percent.  Thus,  the  laminate 
is  subjected  to  a transverse  residual  stress  of  about  10  to  18  MPa. 
Unfortunately,  none  of  the  specimens  was  examined  immediately  after  the 
cure,  and  hence  it  is  not  clear  at  present  whether  or  not  this  stress 
is  sufficient  to  induce  failure  over  a period  of  at  least  a year.  Fur- 
ther work  is  being  done  to  study  the  time-dependent  ply  failure  due  to 
residual  stresses. 

Finally,  the  variation  of  the  maximum  residual  stress  with  the 
equilibrium  moisture  content  is  shown  in  Figure  16  for  various  laminates. 
The  maximum  transverse  stress  is  at  the  mid  surface  for  the  [0^/90^]T 
laminate  and  in  the  0-deg  ply  for  the  [O/OO^l^  laminate.  The  figure 
also  includes  the  longitudinal  shear  and  transverse  stresses  in  the 
45-deg  ply  of  the  [02/±45]j  laminate. 

Complete  removal  of  the  curing  stress  in  all  the  laminates  analyzed 
occurs  at  the  equilibrium  moisture  content  of  approximately  1.12  percent. 
The  corresponding  relative  humidity  can  be  found  by  noting  first  that 
the  equilibrium  moisture  content  is  proportional  to  some  power  of  the 
relative  humidity  H (Reference  1). 

c„  = 

The  parameters  a and  b can  be  calculated  by  using  the  data  in  Reference 
3,  as  follows: 


a * 1.7%  and  b = 1.675 

Thus,  the  relative  humidity  at  which  no  residual  stress  is  present  at 
27°C  is 

H = 78% 
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In  fact,  a simple  relationship  between  the  temperature  and  the 
relative  humidity  for  the  residual  stress-free  state  can  be  obtained 
if  aj  is  neglected.  In  terms  of  the  curing  stress-free  temperature  T^, 
the  relationship  becomes 


With  the  cure  temperature  as  T^,  the  foregoing  equation  is  shown  graphi- 
[ cally  for  the  AS/3501-5  composite  in  Figure  17. 
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SECTION  IV 
CONCLUSIONS 

The  swelling  behavior  of  a graphite/epoxy  (AS/3501-5)  can  be 
described  fairly  well  by  a model  which  is  based  on  the  assumption  that 
the  swelling  is  negligible  until  the  moisture  concentration  reaches  a 
threshhold  value  and  then  increases  linearly  thereafter.  The  parameters 
involved  in  the  experimental  correlation  are  the  transverse  and  thick- 
ness strains  of  symmetric  laminates,  and  the  warpage  of  unsymmetric 
laminate.  The  analysis  also  includes  the  curing  stresses  resulting  from 
fabrication. 

The  void  content  calculated  from  the  threshhold  value  c^  is  0.64 
percent.  It  is  interesting  to  note  that  this  value  falls  in  the  range 
applicable  to  most  graphite/ epoxy  composites  fabricated  under  a standard 
procedure  [11]. 

The  initial  absorption  in  virgin  specimens  seems  to  produce  micro- 
cracks, opening  up  the  passageway  for  the  subsequent  diffusion.  These 
microcracks  are  not  visible  under  x300  magnification  and  entrap  only 
a small  amount  of  additional  moisture. 

The  residual  stresses  appear  to  be  responsible  for  the  absorption 
initially  being  slower  than  the  desorption.  This  is  because  the  boun- 
dary layer  is  in  compression  during  absorption  whereas  it  is  in  tension 
during  desorption. 

Many  cracks  are  found  in  the  unsymmetric  laminates  whereas  none  are 
present  in  the  symmetric  laminates.  At  present  it  is  not  clear  whether 
or  not  this  is  because  of  any  difference  between  AS/3501-5  and  AS/3501-5A. 
In  any  case,  the  processing  parameters  such  as  cool-down  rate,  post  cure 
temperature,  and  temperature  distribution  may  affect  the  state  of  residual 
stress. 
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■ Since  the  present  study  Is  only  a first  step  toward  understanding 

, the  residual  stresses  and  internal  damages,  the  material  behavior  has 

i been  idealized  as  elastic,  and  any  coupling  between  moisture  and  temper 

ature  has  been  neglected.  Therefore,  improved  models  may  be  required 
to  analyze  the  material  response  under  complex  environmental  exposures. 
For  example,  the  effects  of  thermal  spiking  (Reference  12)  may  necessi- 
tate the  development  of  moisture- temperature-dependent  viscoelastic 
models. 
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